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Abstract 

“f?  A unique  application  of  the  Monte  Carlo  method  was  developed  for  determining 
reliability  vs.  cycles  to  failure  of  the  M60  tank  torsion  bar.  In  applying  the 
method,  material  torsional  fatigue  and  spectrum  loads  were  modelled  such  that 
vai lability  in  the  functional  parameters  and  operational  loads  were  represented. 
Random  torsional  displacement  values  obtained  from  the  amplitude  displacement 
distributions  applied  to  the  fatigue  equations  resulted  in  an  exponential  distri- 
bution for  cycles  to  failure  of  the  in  service  bar.  The  number  of  simulations  in 
the  Monte  Carlo  process  was  determined  from  a convergence  criteria  involving 
stability  of  the  third  and  fourth  moments  of  the  cycles  to  failure  distribution. 

Reliability  vs.  bar  life  computations  indicated  a negligible  amount  of  life 
flaw  initiation.  Assuming  a design  change  involving  a twenty  percent 
reduction  in  bar  stresses  increased  the  life  estimates  by  a factor  of  three.  An 
increase  in  reliability  can  also  be  realized  if  computations  are  made  by  assuming 
a bar  has  been  in  operation  for  a specified  number  of  cycleu.  A comparison  of 
minimum  life  (ninety  nine  percent  probability  of  survival)  between  predicted  and 
in  service  results  showed  excellent  agreements  (less  than  eight  percent  difference). 
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Introduction 


The  current  need  for  establishing  reliabiity  of  various  components  and 
systems  for  U.S.  Army  weapon  vehicles  is  being  realized.  The  consequences  of 
over. or  under  design  are  often  reflected  in  either  premature  failure  or  excess- 
ive costs  and  poor  performance  due  to  excessive  weight.  The  mean  life  estimates 
used  as  a criteria  for  defining  acceptability  of  cyclic  loaded  component  will 
often  provide  a false  sense  of  security  regarding  its  capability.  The  applica- 
tion of  higher  strength  ferrous  materials  or  the  less  conventional  structural 
materials  such  as  composites  and  ceramics  will  often  result  in  premature 
failure  because  of  the  inability  to  recognize  the  inherent  variability  of  the 
materials  strength. 


The  objective  of  this  paper  is  to  determine  a methodology  which  will 
circumvent  the  present  deterministic  approach  used  in  establishing  an  acceptable 
design  for  cyclic  random  loaded  structure.  Instead  of  analyzing  the  worst  case 
situation  related  to  the  spectnnn  loads,  S/N  curve,  or  crack  propagation  laws, 
the  authors  introduce  a method  which  simulates  the  variability  in  loading  and 
materials  capability.  Use  of  this  methodology  eliminates  the  over  (worst  case) 
or  under  design  (mean  life)  situation  by  Introducing  a probabilistic  design 
criteria.  Recognition  of  the  reliability  values  as  a function  of  the  life  cycles 
of  operation  can  provide  the  opporttinity  for  selecting  a specified  life  value 
corresponding  to  the  probability  estimate.  The  remaining  component  life  can  then 
be  determined  as  related  to  its  probability  number. 

The  recommended  ASTM  procedure  for  determining  acceptable  design,  involves 
establishing  a lower  confidence  3 Standard  Deviation  bound  on  the  S/N  Curve  then 
selecting  cycles  to  failure  from  the  bounded  curve  consistent  with  predetermined 
maximum  stress  obtained  from  the  spectrum  load  results.  This  procedure  can  often 
result  in  an  over  design  situation  since  the  maximum  load  may  rarely  occur  in- 
addition to  the  fact  there  is  a small  chance  that  the  lower  S/N  Curve  bound  is 
representive  of  the  True  S/N  Curve. 
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The  Monte  Carlo  process  used  in  predicting  life  time  versus  reliability  of  the 
M60  torsion  bars  had  a prior  application  in  a report  by  (1).  Conceptually,  this 
method  is  quite  simple,  requiring  modelling  of  the  spectrum  loads  end  the  material 
fatigue  life  with  respect  to  crack  propagation  or  stress/cycles  to  failure. 

Amplitude  Displacement  Model 

In  figure  1,  a schematic  of  the  torsion  bar  in  the  M60  Tanks  is  shown.  The 
amplitude  distributions  of  three  bars  from  tests  conducted  at  Aberdeen  Proving 
Groiinds  (APG)  is  shown  in  figure  2 • Positive  and  negative  angular  displacements  of 
the  bare  as  function  of  tank  travel  are  shown  in  figure  figure  2b  the  amplitude 

distributions  are  listed  in  a manner  describing  percent  time  less  than  by  a plus 
■ign  (+)  and  percent  time  grater  than  by  a minus  sign  (-) , (eg.  25%  level  equals  a 
-75%  level.  The  + peak  represents  maximum  angular  displacement  under  load,  the 
negative  peak  is  maximum  unloaded  angular  measure.  In  order  to  eliminate 
considering  positive  and  negative  peak  values  in  figure  2a  for  determining  angular 
displacements  in  the  cyclic  loading  process,  the  angular  displacement  is  defined 
as  follows. 
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A0  = e + 1 6" 


where 


8~=  naximun  negative  angular  displacement 


Q _ displacement  from  figure  2b 


d 6 represents  the  adjusted  angular  displacement 


The  Beta  distribution  provided  the  best  representation  of  the  skewed 
^lltude  distribution.  The  dampening  effects  that  occured  under  load  resulting 
from  a stop  used  in  preventing  further  angular  twist  of  the  bar  producing  a highly 
skewed  discrete  cumulative  probability  values.  The  Beta  function  is  defined  as: 


' (r(p)r(Q) 

and  0 s A6  s 1 P,  Q > 0 


The  and  Q values  are  selected  in  a manner  that  provides  the  best  Probabilitv 
Density  Faction  (PDF)  for  representing  the  data.  Figure  3 describes  a typical 
distribution  and  Table  1 shows  the  excellent  correlation  between  predicted 
(Beta  representation)  and  actual  test  result*?  Anoiae  •.  , onO 
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therefore, a good  representation  below  this  angle  is  not  essential. 


Crack  Growth  Law  For  Estimating  Torsion  Bar  Life 
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where  R is  ® uniform  random  number  and  corresponds  to  the  desired  type  of 

frequency  distribution  for  the  parameter.  A probability  density  function  for  the 
N cycles  to  failure  can  be  obtained  by  randomly  selecting  from  C.,  C,,  A , and  A0 
distributions  of  discrete  sets  of  numbers  and  substituting  them  into^equation  3. 
Note«^  there  should  be  an  equal  amount  of  random  numbers  for  each  parfuneter  to 
have  ’proper  amount  of  numbers  for  the  N distribution. 

(B)  S/N  Curve  Analysis 

Torsional  bar  life  expectancy  was  obtained  using  the  Monte  Carlo  process 
applied  to  the  S/N  Curve  relationship.  The  procedure  provided  a method  for 
obtaining  life  time  estimates  of  the  bar  by  combining  the  effects  of  crack 
initiation  and  propagation.  A description  of  the  S/N  Curve  is  shown  in  figure 
6,  where  the  base  line  data  was  obtained  from  a literature i survey  for  material 
metallurgically  similar  to  the  torsion  bar  material.  The  survey  provided  a set 
of  S/N  Curves  for  torsional  fatigue  shown  below  for  hest  representing  the  current 
materials  used  in  the  bar. 

LogloN  * B + .O68A0  (7) 

where  B « 7.70 

The  slope  value  of  .068  was  essentially  the  same  for  all  curve;  in  the  set. 
The  adjustment  in  B from  7.70  to  8.06  made  on  the  basis  of  M60  torsion  bar  quality 
assurance  tests  at  a single  A6  value  performed  at  the  Scranton  manufacturing 
facility  (See  figure  6).  A single  load  equivalent  to  a 42  degree  angular  dis- 
placement was  applied  during  the  quality  assurance  torsional  fatigue  test.  Using 
the  mean  value  and  the  cycles  to  failure  in  Figure  7 provided  a more  accurate 
estimate  of  (B).  The  curves  representing  a range  of  10  and  20  percent  reduction 
in  bar  stress  are  shown  in  figure  6. 

The  S/N  Curve  Monte  Carlo  process  is  similar  to  the  previously  outlined 
method  for  da/dN  relationships.  The  primary  difference  involves  using  Models 
for  (B)  and  A0  from  figure  6 and  2 respectively.  A schematic  of  uhe  basic  S/N 
representation  is  shown  in  figure  8a  and  8b.  In  figure  8a  simulation  of  S/N 
curve  variability  is  shown  for  a specific  value.  Figure  8b  describes  probability 
density  function  (PDF)  for  (B).  A random  selection  of  a discrete  set  of  numbers 
fromA0  and  (B)  distributions  is  then  applied  to  equation  7 in  order  to  obtain 

'^slue.  The  process  is  repeated  until  all  values  from  the  two  distributions 
are  selected.  This  process  will  then  provide  a PDF  to  represent  Log  N. 


Torsion  Bar  System  Reliabilit^ 
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By  assuming  a tank  with  a N torsion  bar  system  the  following  procedures 
would  be  applied  in  order  to  establish  reliability  of  the  system.  If  any  one 
bar  could  cause  failure  (independence)  then  reliability  R will  be 
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Pj  - Prob.  of  Survival 
j - j*"^  Torsion  Bar 
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if  it  is  assumed  that  all  torsion  bars  must  fail  for  system  failure  (dependence) 

then* 


R « ?!  X P2  / Pj  X ...X 


where  ^n-1' 

of  N“1 , 1 bars. 


is  the  reliability  of  Nth  bar,  given  reliabilities 


Reliability  of  Operation  After  Specific  Number  of  Cycles 


The  reliability  of  operating  an  additional  number  of  cycles  when  a 
specified  number  of  cycles  of  operation  has  been  completed  is  obtained  in 
the  following  manner.  Initially  it  is  assumed  that  a specified  distribu- 
tion function  say  f(N)  is  known.  For  example  the  distribution  of  Log^.N 
from  Monte  Carlo  method  previously  described.  The  reliability  R(n  , n) 
is  a conditional  probability  requiring  the  probability  of  operating  for 
+ n cycles  when  n^  cycles  have  been  completed.  That  is 


R(n^  + n) 


R(n^  + n) 


ff(N)dN 


R(n^) 


n^  + n 


*f(N)dN 


where  n is  the  additional  mission  in  cycles  after  n.,  cycles  of  operation. 
The  number  N^  (n^^,  r.)  of  components  (torsion  bars)  feat  will  survive  an 
additional  n cycles  is  given  by 


Ng(n^  n)  » Ng(n^)  . R(n^  n) 


where  N (n.) 
cycles.® 


number  of  components  starting  the  mission  of  n additional 


Results  and  Discussion 


The  proper  number  of  simulations  for  the  Monte  Carlo  Method  depended  on  the 
models  under  consideration.  For  example  5000  and  3000  were  required  for  the 
da/dN  and  S/N  curve  models  respectively.  Using  a convergence  rate  criteria  for 
the  calculated  1 percent  values  (see  Ps  in  figure  9)  and  recognition  of  the  third 


and  fourth  moment  stability  of  the  Log.»N  distribution  provided  an  excellent 
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method  for  determining  required  number  of  simulations.  Differences  in  percentile 
values  for  C.Vc's  of  10  and  15  percent  were  minimum.  The  10  percent  value  was 
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used  for  all  calculations. 
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Although  excellent  agreement  has  been  obtained,  the  authors  would  have 
preferred  representing  the  spectrum  load  consistent  with  an  individual  peak  to 
peak  angular  displacement.  The  simplification  applied  using  the  negative  peak 
as  base  and  representing  the  displacement  relative  to  this  value  was  a g-.'od 
approximation  to  the  available  Individual  displacements.  This  approximation  would 
provide  a slightly  conservative  estimate  in  the  reliability  values r Using  the 
ASTlf-  recommended  practice  of  representing  lower  3 standard  deviation  band  of  the 
S/N  curve  as  measure  of  material  fatigue  loading  capability  combined  with  max- 
imum angular  displacement  (46  degrees)  for  25nph.  The  tank  operation  resulted 
In  a minimum  life  estimate  of  112  miles  for  the  bar.  Selecting  this  number  as  a 
design  allowable  could  result  in  an  overly  conservative  estimate.  The  chance  that 
this  maximum  displacement  could  occur  and  the  S/N  curve  was  the  actual  lower  band 
described  above  is  extremely  small. 

A minimum  life  of  575  miles  was  obtained  from  using  the  maximum  A0  displacement 
value  with  original  S/N  curve  where  .B*  8.06.  This  result  is  obviously  wrong 

limited  samples  of  23  bar  failures  two  of  them  failed  at  mileage  less 
than  400  miles  (See  figure  13).  -i-^eage  less 


Conclusions 


1.  A methodology  for  obtaining  reliability  of  the  M60  tank  torsion  bar 
subjected  to  cyclic  random  loads  has  been  developed  where  probability  of 
survival  is  represented  as  function  miles  of  tank  travel. 

2.  The  developed  methodology  could  be  applied  to  other  structures  with 
cyclic  random  loads. 


3.  The  use  of  the  method  appears  justified  from  recognition  of  the  excellent 
agreement  between  predicted  reliability  estimates  and  those  obtained  from  the 
actual  bar  life  (miles  to  failure)  experienced  during  the  tank  operation. 

4.  Determination  of  minimum  bar  life  was  21  times  greater  from  application  of 
S/N  curve  model  than  that  of  the  assumed  da/dN  model.  This  indicates  most  of 
the  bar  life  exist  prior  to  crack  initiation. 

5.  Application  of  deterministic  procedures,  (use  of  lower  3 5.D  bound  for  S/N 
curve  (ASTM  method)  and  mean  S/N  curve  providing  over  and  under  design  allowable 
estimates  while  Monte  Carlo  method  outlined  in  the  text  values  accurately 
described  acceptable  design  values. 
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Figure  2b 


Detcrminatiaa  of  Qyclea  to  Failure  (ISm  Appr:) 
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gpectrui  Load  (Profile  IV  Course)  - Beta  Function  Representation 


CunaMlative 

Probability 

B (Degrees) 
Test  Results 

B (Degrees) 

Beta 

Representation 

Run  40  (5  nph) 

• 10 

• 14 

.86 

• 25 

4.4 

5.8 

• 14 

5.0 

7.2 

• 50 

8.6 

9.3 

• 66 

12.5 

11.5 

.75 

14.2 

12.7 

.95 

17.0 

16.7 

Run  42  (10  nph) 

.10 

14.0 

6.1 

.25 

16.0 

19.2 

.34 

22.8 

21.5 

.50 

25.8 

24.8 

.66 

29.7 

27.6 

.75 

30.6 

29.0 

.99 

32.6 

32.5 

Run  48  (25  nph) 

.10 

2.3 

10.8 

.25 

22.7 

26.2 

.34 

27.2 

29.1 

.50 

33.7 

33.3 

.66 

39.5 

37.1 

.75 

41.6 

39.3 

.99 

46.0 

46.9 

Cunnulative  Ti»e  Probabilities  of  Tofsional  Bar  Angular  Displacement 
Q adjusted  to  positive  range  by  * 0 where  - «ax. 

negative  angular  displacement. 
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da/dN  Curve  Results 


Mileage  Expected  (Function  of  Spline  Stress) 
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Monte  Carlo  Results  for  S/N  Curve  Mininum  Life 


Estimate 

(99%  Probability  of  Survival) 

VS  Velocity  (MPH) 

Velocity 

Mileage 

Expected 

(MPfi) 

Current  Design 

10%  Design 

20%  Design 

Improvement 

Improvement 

12970 

5 

6,974 

9474 

10 

2,000 

3138 

4420 

15 

345 

638 

1089 

20 

276 

515 

860 

25 

292 

557- 

865 

*Note:  A 99%  survivability  estimate  of  262  miles  was  obtained  from 

cummulative  APG  mileage  on  vehicles  at  time  of  torsion  bar 
failure.  Velocity  of  vehicle  during  tests  was  approximately 
15  to  25  nph. 


Table  3 
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